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(57) Abstract 



A method and apparatus for displaying complex impedance and related physical characteristics of an object (30) at a plu- 
rality of sinusoidal frequencies over a wide frequency range. Excitation signal waveforms (10, 14. 18) are stored in digital form. 
Digital sampling and processing circuitry (84, 86) provide high form. Digital sampling and processing circuitry* (84, 86) provide 
high noise immunity and accurate measurement at all impedance phase angles. The technique is particularly 'well-suited to the 
measurement of complex impedances in living biological tissue because the digital implementation results in highly accurate mea- 
surement values over a plurality of selectable frequencies within a wide sinusoidal frequency region. A novel adaptation of digital 
cross-correlation (84) and convolution (86) techniques is used to simultaneously display (72) real and imaginary electrical imped- 
ance of living biological tissue at several selectable sinusoidal frequencies. 
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METHOD AND APPARATUS FOR DISPLAYING 
MULT I -FREQUENCY BIO-IMPEDANCE 



5 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

My invention relates generally to methods for measurement 
10 and display of complex sinusoidal impedance and, more 
specifically, to the display of complex sinusoidal electrical 
impedance in biological tissues at many frequencies, 

2 • Background of the invention 

15 

The concept of a complex impedance for an object is well- 
known in the acoustical, electrical, and mechanical arts. An 
object impedance is usually defined as the ratio of a 
motivating excitation divided by a resulting object response 

2 0 at a single sinusoidal frequency, gj. For example, the complex 
mechanical, impedance, Z, of a structure can be expressed as 
the ratio of a sinusoidal excitation force, F, to the 
resulting sinusoidal velocity, V, of the structure at the 
point of application of the excitation force, that is: 

25 2(w) = F(w,t) /V(w, t) 

The complex nature of such a driving point impedance 
arises from the time delay, d, of the peak sinusoidal 
velocity, V, with respect to the peak sinusoidal excitation 
force, F. This is commonly expressed in the following form: 
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r - r • i-Y- ^ — 2 

or ~ = -r- (cos cod + i sincod) = 



Where: 

a: « angular sinusoidal frequency in radians/second, and 

cad = phase delay of V with respect to F in radians. 

This concept of complex object impedance at a single 
sinusoidal frequency can be used to express a mechanical ratio 
of force to velocity, an acoustic ratio of pressure to 
displacement, an electrical ratio of voltage to current, a 
thermal ratio of temperature differential to heat flow, an 
electromagnetic ratio of electric field to magnetic field, and 
so forth as is well-known in the art. 

The classical method known in the art for measuring and 
displaying values of complex object impedance requires the 
measurement of the response of the object to an applied 
excitation at a single sinusoidal frequency. For instance, 
the complex electrical impedance of an object can be 
determined, by applying a sinusoidal voltage to the object and 
measuring the resultant sinusoidal current flow through the 
object. The electrical object impedance magnitude may then be 
determined as the ratio of the root-mean-square (RMS) voltage 
and current values and the object impedance phase angle may be 
determined as the delay in radians of the peak sinusoidal 
current with respect to the peak sinusoidal voltage. The real 
and imaginary components of the complex object impedance may 
be determined from impedance magnitude and phase ancle 
according to: 
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Resistance = Re[Z] = R = |z|cos£, where 6 = cod, 
Reactance = Im[Z] * X = |Z|sin0, and 
Impedance = Z = R + iX 
A number of problems are known in the impedance 
5 measurement art that lead to errors and ambiguities when the 
values of complex object impedance are measured and displayed. 
The most significant problem arises from the requirement that 
the complex impedance be determined at a single sinusoidal 
frequency. The presence of other frequency components in the 
ro excitation and response signals lead to errors in the 
displayed complex impedance. Practitioners in the art have 
proposed analog filtering techniques and special excitation 
signal generation methods for minimizing the errors arising 
from such unwanted signal components. However, analog filter 
15 devices are subject to calibration errors, thermal drift and 
other problems, which introduce amplitude and phase errors in 
the excitation and response signals. These phase errors often 
become severe at high values of reactance (imaginary 
impedance) and can completely overwhelm the character of the 

2 0 complex impedance under measurement. 

Other serious limitations to accurate complex impedance 
display are well-known in the art. High impedance magnitudes 
require the determination of a ratio between a large 
excitation signal and a very small response signal, leading to 
25 errors arising from the presence of noise in the small 
response signal . Practitioners propose the use of very high 
excitation signal amplitudes to overcome the effects of such 
noise, but high levels of excitation signal may damage or 
destroy the object to which the excitation signal is applied. 

3 0 This is especially problematic in cases where the objecr under 

measurement is living biological tissue. 
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Another problem arising from the single frequency nature 
of complex object impedance is the difficulty in measuring 
complex cb jeer impedance at several frequencies over a wide 
frequency range. Because signal filters are necessary to 
ensure sinusoidal purity, these filters must be retuned to 
permit impedance display at other sinusoidal frequencies. The 
range of frequencies over which analog filters can be retunec 
is very limited , often to a few octaves. Alternatively, 
individual signal filters can be provided for each measurement 
frequency, but this approach seriously limits the number cf 
frequencies for which complex object impedance can be 
displayed in any particular situation. 

The analog signal filtering and display components used 
in the art are subject to calibration drift resulting from 
thermal changes and variations in operating region. This 
problem leads to a requirement for frequent recalibration to 
minimize compl.ex impedance display errors, thereby preventing 
the rapid and effective measurement of complex object 
impedance at several frequencies. 

The use of balanced bridge impedance measurement 
techniques known in the art can overcome many cf the 
deficiencies of the signal ratio measurement methods mentioned 
above. However, balanced bridge techniques require a well 
calibrated impedance standard and are limited in practical 
application to the determination and display of electrical 
impedances. Moreover, the standard analog bridge components 
are presumed to be linear with respect to signal amplitude, 
which is an inaccurate presumption in most applications. This 
presumption leads to errors in display of complex object 
impedance that arise with changes in excitation signal level. 
Finally, balanced bridge techniques are unsuitable for 
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accurate complex impedance determination over a wide frequency 
range and, in fact, are often limited to a single sinusoidal 
calibration frequency at a single excitation signal amplitude. 
Another problem well-known in the art arises from the 
5 effects of small errors in impedance phase angle at angle 
values approaching n/2 radians (90°) . At these phase angles, 
the tangent function approaches infinity and extremely small 
errors in phase angle cause extremely high errors in one or 
both of the complex impedance components. A number cf clever 

10 techniques have been proposed by practitioners in the art tc 
overcome this problem, but most rely on the thermal 
discrimination between real and imaginary electrical power 
flow in a circuit and are unsuitable for use over a range of 
frequencies or with nonelectrical impedance measurement . 

15 Yet another problem with classical impedance measurement 

and display techniques is the difficulty in simultaneously 
determining impedance at widely separated frequencies. This 
problem arises when the object under measurement experiences 
changes in impedance properties as a result of the application 

2 0 of the excitation signal. An example of this is the well- 
known tendency of an object to increase in temperature in 
response to the application of an electrical voltage. When 
the complex object impedance is a function of object 
temperature , then the measurement of complex impedance at one 

25 frequency will heat the object and create errors in the 
measurement of complex impedance at another frequency because 
of the difficulty in measuring such impedances simultaneously • 
One solution is to use a plurality of impedance measurement 
and display devices to simultaneously measure complex object 

20 impedance at a plurality of sinusoidal frequencies, but this 
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nethod is cumbersome and not practical for large numbers cf 
measurement frequencies. 

These and other related difficulties with accurate 
display of complex object impedance are exacerbated in the 
situation where the object of interest is living biological 
tissue. The measurement and display of complex electrical 
impedance in biological tissue has been of interest since the 
late 1800' s. By 1921 it had been well-supported that the 
living cell had a well-conducting interior surrounded by a 
relatively impermeable, poorly conducting membrane. In 1S25, 
Fricke reported (Fricke, H. , Mathematical Treatment of the 
Electrical Conductivity a nd Capacity of Disperse Systems, 
Phys. Rev. 26:678-681, 1925) that at low frequency (LF) 
currents, there was little conduction through the cell because 
of high membrane capacitance. Thus, Fricke argued that 
conduction occurs primarily in the extracellular fluid (ECF) 
compartment and, at a sufficiently high frequency (HF) , the 
current is shunted through the cell membrane and conducts 
through both the ECF and intracellular fluid (ICF) . 

Knowledge of biological material properties has been 
obtained through complex impedance measurements across various 
cells, suspensions, fibers, eggs and tissues. Of primary 
significance were the discoveries of cell membrane 
capacitance, the beta dispersion region, the additional 
regions of dispersions for cell suspensions, and Maxwell's 
mixture theory for analyzing the impedance measurement data. 

In general, all cells and tissues may be expected to show 
three major dispersion regions in relation to frequencies (a, 
3 and 7) . Of particular interest is the central ^-dispersion 
region, which is explained by the dielectric capacity cf the 
cell membranes. The a and 7 regions are attributed to a 
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surface conductance and to intracellular components. The 3 
region can be expressed in either a simple or complex 
equivalent circuit orientation form, where parallel resistance 
(R P ) reflects extracellular fluid (R ww ) and the series 
5 orientation corresponds to intracellular fluid resistance 
(R icw ) and cell membrane capacitance (C M ) . The equivalent 
circuit used in the Maxwell analysis must give the semi- 
circular relation (impedance and admittance loci) between the 
real and imaginary components of the impedance as applied 

10 frequency is varied, according to Cole (Cole, K.S. , Membrane , 
Ions and Impulses/ University Press, Berkeley, 1968) . . - 

An early attempt to apply this concept was made by 
Thomassett (Thomassett, A. , Bio-electrical Properties of 
Tissues . Lyon Med. 209:1325-1352, 1963), who measured simple 

15 impedance at high (100kHz) and low (1kHz) frequency currents 
using a two-wire technique- The high correlations of high 
frequency impedance to total body water (TBW) and low 
frequency impedance to extracellular fluid gave further 
support to the membrane and dispersion theories- However, 

20 Thomassett 1 s approach did not develop into a practical 
clinical method. 

Using a four-wire configuration, Hoffer, et al. (Hoffer, 
E. / Meador, C. , Simpson, C. , Correlation of Whole-Bodv 
Impedance with Total Body Water Volume . J. Appl. Physio. 17, 

25 4:531-534, 1969) reported a high correlation between measured 
TBW and TBW estimated by Ht 2 /Resistance at lOOJcHz. However, 
Hoffer, et al. reported that their high standard deviations 
implied that further development was necessary to make the 
technique a practical clinical method. 

30 A decade later, Nyboer, et al. (Nyboer, J., Liedtke, 

Reid, K. , Gesert, w., Nontraumatic Electrical Detection of 
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Total Body Water and Density in Man . Proceedings of the VI 
ICEBI, 381-384, 1983) found that measurements of electrical 
resistance at 50kHz, combined with subject weight, height and 
age, could accurately determine body density (fat and lean) in 
human subjects. These relationships were based on the 
assumption that a strong relationship must exist between Fat- 
Free Mass (FFM) and TBW estimated by impedance because FFM 
tissue is consistently hydrated. Following the Nyboer, et 
ai. work, a tremendous amount of interest arose for the 
impedance method. 

Since the Nyboer, et al. presentation, many practitioners 
have sought to validate the single frequency impedance method 
of estimating human body composition in various populations 
(Lukaski, H.C., Johnson, P.E., Bolonchuk, W.W. , Lykken, G.I., 
Assessment of Fat-Free Mass Using Bioelectrical Impedance 
Measureme nts of the Human Body, Am. J. Clin. Nutr. , 41:810- 
817, 1985; Segal, K. , Van Loan, M. , Fitzgerald, P., Hodgdon, 

J. A., Van Itallie, T.B., Lean Body Mass Estimation by 

Bioelectrical Impedance Analysis: a Four Site Cross- 
validation Study , Am. J. Clin. Nutr. 47:7-14, 1988; and 
Kushner, R. , Kunigk, A. , Alspaugh, M. , Andronis, F.,Leitch, 
C. , Schoeller, D. , Validation of Bioelectrical-Inoedance 
Analysis as a Measurement of Change in Body Composition in 
Obesity , Am. J. Clin. Nutr. 52:219*23, 1990). Results have 
been mostly positive for normal subjects but the predictions 
of FFM have been far less precise in the clinical and abnormal 
population groups. Although impedance data continued tc 
correlate highly to TBW, the underlying assumption that FFM is 
consistently hydrated has recently been found to be incorrect 
(Deurenberg, P. , Westtrate, J. A., Hautvast, J., Changes in 
Fat-Free Mass During Weight Loss Measured by Bioelectrical 
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Impedance and bv Densitometry , Am. J. Clin. Nutr. 49:33-6, 
1989) . Furthermore, increased awareness of the usefulness of 
the multi-frequency measurements has led to a general belief 
that the single frequency method is too simplistic and 
5 limiting (Cohn, S., Hov Valid are Bioelectr ical Impedance 
Measurements in Body Composition Studies? . Am. J. Clin. Nutr. 
42:889-890, 1985). Criticism was also raised over the 
clinical usefulness of TBW in view of frequently encountered 
fluid shifts between the ICW and ECW. Van Itallie, et al. 

10 (Van Itallie, T., Segal, K. , Nutritional Assessment of 

Hospital Patients: New Methods and New Opportunities , Am. J. 
Hum. Bio 1:205-8, 1989) assert that a practical method of 
discerning the ICW from the extracellular fluid would offer 
much greater utility and could profoundly influence hospital 
15 patient care and diagnosis. 

Recently, several practitioners (Lukaski, H.C., 
Bolonchuck, W.W. , Estimation of Body Fluid Volumes Using 
Tetrapolar Bioelectrical Impedance Measurements , Aviat. Space 
Environ. Med. Dec, 1163-69, 1988 and McDougall, D. , Shizgal, 
20 H., Body Composition Measurements from Whole Body Resistance 
and Reactance , surgical Forum, 36:43-44, 1986) have proposed 
using the -reactive element in a complex single frequency 
(50kHz) impedance measurement to accurately discriminate the 
extracellular from the cellular mass. However, the use of 
25 multi-frequency measurements of impedance remains the 
technique of choice (Boulier, A., Fricker, J., Thomassett, 
A.L., Apfelbaum, M. , Fat-Free Mass Estimation by Two-Electrode 
Impedance Method , Am. J. Clin. Nutr. 52:581-5, 1990). 

Several practitioners have continued to test the 
30 Thomassett technique (Jenin, P., Lenoir, J., Roullet, C. # 
Thomassett, A., Ducrot, H. , Determination of Body Fluid 
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Compartments by Electrical Impedance Measurements , Aviat . 
Space Environ. Med. 46:152-5, 1975; Settle, R.G., Foster, 
K.R., Epstein, B.R. , Mullen, J.L., Nutritional Assessment : 
Whole Body Impedance and Body Fluid Compartments , Nutr. 
Cancer, 2:72-80, 1980; and Tedner, B.T., Equipment Using an 
Impedance Technique for Automatic Recording of Fluid-Volume 
Changes During Hemodialysis , Med. & Biol. Eng. & Comput. 
21:285-290, 1983) . These practitioners use the ratio of high 
to low frequency simple impedance to reflect fluid compartment 
volume, the normal and abnormal fluid ratio, and fluid 
compartment change. The major problem reported for this 
technique is that the simple HF/LF ratio is too simplistic 
because volume is not determined and the compartment actually 
affected is not identified, although a change in ratio does 
reflect a change in compartmental volume. Furthermore, 
practitioners note that the lack of two-frequency simple 
impedance instrumentation prevents further exploration of this 
approach . 

Because living tissue is mainly affected by ECF at 
frequencies below the ^-dispersion region, and ECF, C M and ICF 
at frequencies above this region, Kanai, et al. (Kanai, K. , 
Haeno, M. f Sakamoto, K. , Electrical Measurements of Fluid 
Distribution of Legs and Arms , Med. Prog. Tech. 12:159-170, 
1987) computed the component values of the human muscle tissue 
equivalent circuit model (R^, R^, and C M ) known in the art. 
By mathematically analyzing the complex impedance measurements 
at multiple frequencies, Kanai, et al. obtained information 
specific to the ECW and ICW. However, others have been unable 
to replicate this advanced approach because of the lack of 
necessary and appropriate complex impedance instrumentation . 
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Expanding the prior art human body composition model from 
a simple two-compartment (lean/ fat) model to a more complex 
model including ICF and ECF, protein, mineral and fat, creates 
a clearly-felt need for new and improved assessment 
techniques. This need has been unmet, until now because of 
the lack of appropriate complex impedance instrumentation and 
the lack of the necessary degree of sophistication in the 
associated measurement and analysis methods. This situation 
has inhibited the progress of this promising technology. 

Increased use of diuretic drugs, fluid monitoring 
difficulties in intensive care, and the shrinking cell mass 
and expanding ECF that accompany most systemic wasting 
diseases and malnutrition has led to a well-recognized' need 
for effective multi-frequency bioimpedance instrumentation in 
support of the compartmental approach to body fluid 
assessment. These unresolved problems and deficiencies are 
clearly felt in the art and are solved by my invention in the 
manner described below. 

SUMMARY OF THE INVENTION 

My invention is a method and apparatus for the 
simultaneous determination and display of complex electrical 
bio-impedance at several frequencies. Although I teach 
primarily the application of my invention to the display of 
multi-frequency electrical bio-impedance, my invention is also 
readily applicable to the determination and display of complex 
impedances for many other types of objects, including other 
impedances such as mechanical, acoustic and the like. This 
can be appreciated by recognizing that my invention determines 
and displays complex object impedance by treating an 
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excitation signal and a response signal as two independent 
functions in the time domain. Because of this view, Z car. 
apply certain signal processing techniques to the two 
independent time domain signals to obtain results in a manner 
previously unsuspected in the art. 

The first such signal processing concept that I use is 
the creation of a complex cross-correlation signal as a 
function of time delay, t, between the excitation, e, and 
response, r, signals. My complex cross-correlation function . 
r„(t,, can be developed from the integral over time of the 
product of the two excitation and response signals where the 
response signal is delayed by time r with respect to the 

excitation signal. 

The second signal processing concept that I use is the 
complex Fourier transform, F, of the cross-correlation signal, 
r„(t). The Fourier transform of ^(t) is F (R^) , which is a 
signal in the frequency domain having a complex value at each 
angular frequency u that is proportional to the complex object 
impedance . 

Following the conversion of the two excitation and 
response signals to an intermediate cross-correlation signal 
and a final Fourier transform signal, I display the complex 
object impedance as a function of frequency over a 
predetermined range. My cross-correlation signal and Fourier 
transform signals are both complex because they exhibit 
amplitude and phase information at each point in the time 
delay and frequency domains, respectively. 

My application of these simple signal processing 
techniques fcr the conversion of time domain signals tc 
frequency domain signals can be embodied ir. a an apparatus 
using either analog or digital electronic components. Ar. 
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advantage of my preferred embodiment is that my use of digital 
components leads to improved storage and processing accuracy 
and minimizes the effects of analog component calibration 
drift. I have also discovered that the development of a 
5 complete cross-correlation signal, R(t), over a large time- 
delay region and the subsequent development of a Fourier 
transform signal, F(o>) , over a broad frequency region, while 
possible, is not the most useful embodiment of my invention. 
Accordingly, I have refined this elementary convolutior. 

10 approach by adding the following improvements. 

First, I use an electrical current for the excitation 
^signal instead of the electrical voltage normally used in the 
art. The constant current excitation signal prevents 
unsuspected hazards to iiving biological tissues because the 

15 biological hazards are more closely related. to current levels 
than to voltage levels as is known in the medical art. 
^Secondly, I make use of a series of stored digital data to 
define the excitation current waveform precisely. By 
precisely defining the excitation current waveform, I can 

2 0 create and apply a single-frequency excitation current to the 
living biological object under measurement. Conversion of two 
such single-frequency excitation and response signals to a 
cross-correlation signal results in a sinusoidal cross- 
correlation signal having an amplitude and phase 

25 representative of the complex cross-correlation between the 
excitation and response signal. Moreover, the Fourier 
transform of such a sinusoidal cross-correlation signal is a 
single impulse function at the point in the frequency domain 
equivalent to the single frequency of the initial excitation 

3C and response signal. 
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These observations are particularly useful because they 
result in simplification of my signal and convolution 
conversion process to a simple signal multiplication procedure 
in quadrature. That is, the complex object impedance at a 
single frequency can be developed from the excitation and 
response signals in the time domain by a multiplication and 
integration process in quadrature. This quadrature signal 
conversion method is the preferred embodiment of my invention 
because of the relatively simple apparatus necessary. 

A potential problem with my quadrature technique is thar 
a different excitation waveform is required for each 
sinusoidal frequency value at which complex object impedance 
must be displayed. I have resolved this problem by providing 
for the storage of a plurality of excitation waveforms in a 
signal generator digital memory means so that the excitation 
signal can be quickly stepped through a variety of sinusoidal 
frequencies merely by selecting the necessary waveform data 
from the memory device in the desired sequence. The 
excitation current is then produced from the stored digital 
data and the analog response voltage developed across the 
object under measurement is sensed. Moreover, network 
analyzer means may be added to compute and display a plurality 
of equivalent circuit elements from a similar plurality of 
complex impedance data obtained at the different sinusoidal 
frequencies stored in my digital signal generator memory 
means . 

For ease of implementation and minimization cf 
measurement and display errors, I prefer to perform all signal 
conversions using digital means. Accordingly, in my preferred 
embodiment, tne response and excitation signals are converted 
tc digital form by sampling means and the multiplication cf 
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the two excitation response signals in the time domain is 
performed using digital multiplication means. The real or 
resistive component of the complex object impedance is related 
to a first impedance signal expressed as the arithmetic 
5 average of the products of a series of simultaneous samples of 
the excitation signal and response signal. The imaginary or 
reactive component of the complex object impedance is related 
to a second impedance signal, expressed as the arithmetic 
average of a series of products of simultaneous samples of the 

10 response signal and a delayed excitation signal, where the 
delay is precisely equal to n/2 radians or one-fourth cycle of 
the sinusoidal waveform representing the frequency at which 
.the complex object impedance is determined. 

By using digital sampling techniques and synchronizing 

15 the sampling windows of the excitation and response signals, 
the implementation of my preferred embodiment is simplified 
and the signal conversion process made mpre robust (resistant 
to noise and calibration errors) . It is not necessary to 
independently sample the excitation current waveform because 

2 0 this waveform is developed from digital samples stored in a 

digital memory device and these samples can be used directly 
in the quadrature multiplication process. However, I prefer 
to independently sample the analog excitation current waveform 
at the point where it is applied to the biological object for 
25 two reasons. First, this simplifies the synchronization of 
the sampling windows for the two excitation and response 
signals. Secondly, this avoids errors arising from unexpected 
changes in the excitation signal between the memory storage 
means and the point of application to the object under 

3 0 measurement. 
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One cf the advantages of my invention is that complex 
bio-impedance can be measured almost simultaneously at a 
plurality of frequencies in a wide frequency region. Another 
advantage of my invention is that a change in the sinusoidal 
excitation signal frequency can be achieved without changing 
analog circuit paths or component values. This minimizes the 
errors arising from component calibration errors. An 
objective of my invention is to avoid complex impedance errors 
arising from nonlinearities in the object under measurement. 

Yet another advantage of my invention is the rejection cf 
noise effects over a broad frequency range without the use of 
analog filter techniques. It is another object of my 
invention to accurately display the complex object impedance 
phase at all phase values from -90° to +90°. Another 
advantage of my invention is that a variety of calibration 
factors may be readily applied to the displayed results by 
using digital techniques known in the art. 

An important objective of my invention is to permit the 
measurement of bio- impedance at a plurality of frequencies in 
as little as two cycles of the primary frequency waveform. 
Another important advantage of my invention is that my time 
series convolution process eliminates most errors that 
otherwise arise from the undesired presence of random noise in 
rhe relevant signals. 

The foregoing , together with other features and 
advantages of my invention will become more apparent when 
referring to the following specifications, claims and the 
accompanying drawings . 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a more complete understanding of my invention, 
reference is now made to the following detailed description of 
5 the embodiments illustrated in the accompanying drawings 
wherein: 

Figure 1, comprising Figures 1(a) through 1(h) , shows a 
series of waveforms in the time domain exemplifying the 
excitation signal, delayed excitation signal, response signal, 
10 first and second impedance signals and the related series of 
digital samples; 

Figure 2, comprising Figures 2A and 2B, illustrates a 
method for applying an excitation signal to an object and 
detecting a resulting response signal and shows a typical 
15 equivalent circuit representation of a living cell; 

Figure 3 is an illustrated embodiment of an apparatus for 
measuring and displaying the complex impedance of an object ; 

Figure 4 is an illustrated embodiment of the signal 
generating means from Figure 3; and 
2 0 Figure 5 is an illustrated embodiment of an apparatus 

that uses my more general convolution method to display the 
complex impedance of an object over a region in the frequency 
domain, 

2 5 DESCRIPTION OF THE PREFERRED EMBODIMENT 

Figure 1 illustrates the typical waveform characteristics 
for the excitation and response signals used in my invention 
for displaying complex object impedance. Figure 1A shows an 

3 0 excitation signal current 10 that varies in amplitude 

sinusoidally with time in the manner shown. Figure IB shows 
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rhe resulting series of excitation pulses or samples 12 thar 
represenr current 10 samples made at a predetermined sampling 
rate. This predetermined sampling rate should be high enough 
(as seen from the Nyquist criteria, known in the art) to avoid 
losing significant phase and amplitude information contained 
in current 10. Figure 1C illustrates a delayed excitation 
signal current 14 that is identical to current 10 except for 
a time delay of rr/2 radians or one-quarter cycle of the 
sinusoidal waveform. Figure ID illustrates a series cr 
delayed excitation pulses or samples 16 that represent current, 
waveform 14 samples taken in synchronism with excitation 
current samples 12 . 

In operation, my invention uses excitation signal current: 
10 to excite a living biological tissue sample and I then 
monitor the sample for a response signal voltage waveform 18 
shown in Figure IE. Response signal voltage IS may contain 
several sinusoidal frequency components, but the fundamental 
sinusoidal frequency of signal current 10 will predominate. 
Figure IF illustrates a series of response voltage samples 2C 
that are taken in synchronism with samples 12 and 16. 

Figures 1G and 1H illustrate the first impedance samples 
22 and second impedance samples 24, respectively. Sample 22 
is obtained by multiplying the concurrent values of sample 12 
and sample 16 (I n V B ) , either by analog or digital means. 
Sample 24 is obtained by multiplying concurrent values of 
delayed excitation sample 16 and response voltage sample 2C 
(l90n v n) • Figure 1G shows a first mean impedance signal 26 
representing the long term arithmetic average cf impedance 
samples 22. Similarly, Figure 1H shows a second mean 
impedance signal 28 representing the long tern arithmetic 
average cf samples 24. 
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While I believe that my invention is useful for 
determining many different types of impedances, including 
thermal, mechanical, acoustical and so forth, I limit my 
discussion herein to my preferred application, which is the 
5 measuring of complex electrical impedances of living 
biological tissues* Figure 2A illustrates the fundamental 
arrangement required to measure and display complex electrical 
impedance of an object 30, Excitation signal current 10 is 
applied to object 30 and response signal voltage 18 is 

10 monitored- 

Figure 2B illustrates a useful equivalent circuit model 
for the living human cell 32 proposed by Kanai, et al. , cited 
above. Cell model 32 comprises three lumped constant RC 
circuits 34, 36 and 38. RC circuit 34 represents the cell 

15 membrane capacitance and cell membrane resistance, RC circuit 
3 6 represents the intracellular fluid capacitance and 
resistance. RC circuit 38 represents the extracellular fluid 
capacitance and resistance. 

Human cell equivalent circuit 3 2 presents a single 

20 complex impedance to the pair of terminals 40 at any single 
sinusoidal frequency. Accordingly, the six lumped-constant 
model elements can be determined only from impedance data at 
six (or more) different sinusoidal frequencies by means of a 
critically-or over-determined system of linear equations as 

25 known in the art. 

As discussed above, I prefer to use a simplified signal 
multiplication procedure in quadrature to determine complex 
object impedance. This signal multiplication procedure was 
introduced and described in connection with Figure 1, wherein 

30 first and second mean impedance signals 26 and 28 represent 
real and imaginary components (Re[ Z] , Im[Z] ) . Figure 3 shows 
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a functional implementation of my invention, including the 
necessary means for determining a plurality of equivalent 
circuit element values such as described in Figure 2E fcr 
human cell 32. These equivalent circuit analyzer elements are 
5 shown connected by dotted lines in Figure 3, and are further 
described below. 

A frequency selection means 42 may comprise front panel 
switch controls, digital memory means or other related devices 
suitable for selecting one or more sinusoidal operating 

1C frequencies. Each selected frequency is passed to a signal 
generator means 44 and also to an impedance analyzer memory 
means 46. Impedance analyzer memory means 46 merely receives 
the frequency selection and stores it in memory for later use 
with concurrent complex body impedance data. Signal generator 

15 means 44 preferably comprises a frequency sequencing means 48/ 
a frequency memory means 50 and a generator memory means 52 as 
shown in Figure 4. Frequency sequencing means 48 selects from 
among a plurality of operating frequencies stored in memory 
means 50 and 52. Frequency memory means 50 provides the 

ZZ necessary clocking pulses to generator memory means 52, which 
then provides a train of excitation current pulses making up 
a sinusoidal excitation current in response to the commands 
received from frequency sequencing means 48 and frequency 
memory means 50. 

25 Returning to Figure 3, a current limiting means 54 is 

interposed in series with excitation current samples 12 to 
prevent the excitation current amplitude from exceeding 
biologically safe levels. Excitation current samples 12 are 
then imposed on object 3 0 at terminals 4 0 and a response 

2Z detector means 56 is connected across terminals 40 to measure 
response signal voltage 18. 
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First and second analog-to-digital (A/D) conversion means 
58 and 60 serve to convert excitation signal current 10 and 
response signal voltage 18 from analog form to digital form in 
a manner known in the art. The operation of A/D means 58 and 
5 60 is synchronized with the operation of signal generator 
means 44 by clocking means (not shown) to ensure that the 
leading and lagging edges of the digital pulses created by the 
three circuits are synchronous in time to within less than one 
microsecond. 

10 The output of first A/D means 58 presents excitation 

samples 12 to the input of a quadrature delay means 62 and the 
first input of a first multiplier means 64. The output of 
. second A/D means 60 presents response samples 2 0 to the second 
input of first multiplier means 64 and the first input of a 

15 second multiplier means 66. The output of quadrature delay 
means 62 provides delayed excitation samples 16 that represent 
a digital embodiment of delayed excitation signal current 14 
discussed in connection with Figure 1. Delayed excitation 
signal current 14 lags excitation signal current 10 by 

20 precisely one-fourth of a waveform cycle at the fundamental 
sinusoidal operating frequency. Delayed excitation samples 16 
from quadrature delay means 62 are presented to the second 
input of second multiplier means 66. 

The output of first multiplier means 64 provides a series 

25 of pulses equal to the series of products of concurrent 
excitation samples 12 and response samples 20, which we have 
named first impedance samples 22 in Figure 1. Similarly, the 
output of second multiplier means 66 provides a series of 
pulses named second impedance samples 24 , which is the series 
3 0 of products of concurrent delayed excitation samples 16 and 
response samples 20. First impedance samples 22 are presented 
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to a first: averaging means 68 and second impedance samples 24 
are presented to a second averaging means 70. The outputs 
from first and second averaging means 68 and 7 0 comprise firsr 
and second mean impedance signals 2 6 and 28 as discussed above 
5 in connection with Figure l. Signals 26 and 28 are exactly 
proportional to the real and imaginary components of the 
complex electrical impedance of object 30. Signals 2 6 and 28 
are presented to impedance display means 72, which is 
preferably a multi-digit LCD numerical display cr ether 

1C suitable display means known in the art. 

The calibration means 74 comprises circuitry that 
produces a current calibration signal 76 and a voltage 
calibration signal 78 for use in adjusting impedance display 
means 72 in response to independent calibration of excitation 

15 signal current 10 and response signal voltage 18. 

Impedance display means 72 also presents the complex 
impedance of object 3 0 to impedance analyzer memory means 3 6 
for storage together with the relevant sinusoidal frequency 
value presented by signal frequency selection means 42. These 

20 data are stored in impedance analyzer memory means 46 for use 
in computing a plurality of lumped constant equivalent circuit 
elements -as discussed above in connection with Figure 2. 
These equivalent circuit elements are computed and displayed 
by means of a network analyzer means 80 and an equivalent 

2 5 circuit display means 82 disposed substantially as shown in 
Figure 3. Network analyzer means 80 acts to solve a system of 
linear equations for a plurality of equivalent circuit 
elements from a plurality of complex impedance and operating 
frequency data pairs in any suitable digital or analog manner 

2: known in the art for solving linear systems of equations. The 
resulting equivalent circuit element values are then presented 
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by network analyzer means 8 0 to equivalent circuit display 
means 82 for display in any suitable manner known in the art 
for displaying data. 

Although I prefer the simple digital quadrature 
5 multiplication and integration process discussed above, I also 
show in Figure 5 the use of a cross-correlator means 84 and a 
convolver means 86 accomplish a similar result in lieu of the 
apparatus disclosed in Figure 3. The effect of cross- 
correlation of signals 18 and 10 followed by convolution cf 
10 the cross-correlation function R*(t) shown in Figure 5 was 
briefly discussed above and is further treated in detail by K, 

G. Beauschamp and C. Yuen, Digital Methods for Signal 

Analysis , George Allen & Unwin, Ltd, London, 1979. " The 
applications of these digital signal processing techniques to 
15 the measurement of complex body impedance is suggested herein 
for the first time. Note that the two outputs 126 and 12 8 
from convolver 86 contain complex impedance data over the 
entire frequency domain and impedance display means 7 2 can be 
Configured to display impedance data at a single frequency in 
20 response to an input from frequency selection means 42. 

Obviously , other embodiments and modifications of my 
invention will occur readily to those of ordinary skill in the 
art in view of these teachings. Therefore, my invention is to 
be limited only by the following claims which include all such 
25 obvious embodiments and modifications when viewed in 
conjunction with the above specification and accompanying 
drawings . 

I CLAIM: 
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CLAIMS 

1. A method for displaying the complex impedance of an 
2 object ar a sinusoidal frequency comprising the steps of: 

selecting a sinusoidal frequency value; 
4 selecting a first sampling rate no less than 4 00% of said 

s inuso ida 1 frequency ; 
6 selecting a second sampling rate no less than first said 

sampling rate; 

c generating a series of excitation pulses at said first 

sampling rate to form an excitation signal comprising a 
10 component at said sinusoidal frequency having a predetermined 
magnitude; 

12 applying said excitation signal to said object; 

detecting a response signal comprising the response of 
14 said object to said applied excitation signal; 

sampling said excitation signal at said second sampling 
16 rate to form a series of excitation samples ; 

sampling said response signal at said second sampling 
18 rate in synchronism with said excitational sampling step to 

form a series of response samples; 
2 0 delaying said excitation samples by a delay time 

substantially equal to one-fourth cycle of said sinusoidal 

22 frequency to form a series of delayed excitation samples; 

multiplying each said excitation sample by the 
2 4 simultaneous said response sample to form a first impedance 
sample ; 

2 5 multiplying each said delayed excitation sample by the 

simultaneous said response sample to form a second impedance 

23 sample; 
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averaging said first and second impedance samples over 
3 0 time to form first and second mean impedance signals, 

respectively; and 
3 2 displaying said first and second mean impedance signals 

to represent said complex object impedance. 

2. The method described in Claim 1 wherein: 

2 said complex impedance comprises the electrical impedance 

of said object; and 
4 said sinusoidal excitation signal comprises an electrical 

current of substantially invariant sinusoidal magnitude. 

3. The method described in Claim 2 wherein: 

2 said first sampling rate is equal to said second sampling 

rate and said first and second sampling steps are synchronized 
4 with said generating step. 

4. The method described in Claim 2 wherein said object 
2 comprises living biological tissue. 

5. A method for displaying the complex impedance of an 
2 object over a range of sinusoidal frequencies comprising the 

steps of: 

4 selecting a first sampling of no less than 4 00% of the 

maximum said sinusoidal frequency; 
6 selecting a second sampling rate no less than said first 

sampling rate; 

8 generating a series of excitation pulses at said first 

sampling rate to form an excitation signal; 
10 applying said excitation signal to the object; 
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detecting a response signal comprising the response of 
said object to said applied excitation signal; 

sampling said excitation signal and said response signal 
in synchronism at said second sampling rate to form a series 
of synchronous excitation samples and response samples, 
respectively ; 

computing the complex cross-correlation of said series of 
excitation samples and said series of response samples over a 
predetermined time delay region; 

computing the real and imaginary components of the 
Fourier transform of said complex cross-correlation over said 
range of sinusoidal frequencies to form first and second 
impedance signals, respectively; 

computing the square of the magnitude of the Fourier 
transform of said excitation signal over said range of 
sinusoidal frequencies to £orm a calibration signal; 

displaying the ratios of said first and second impedance 
signals to said calibration signal to represent said complex 
object impedance over said range of sinusoidal frequency. 

6. The method described in Claim 5 wherein: 

said first sampling rate is equal to said second sampling 
rate and said first and second sampling steps are synchronized 
with said generating step. 

7. The method described in Claim 5 wherein: 
said object comprises living biological tissue. 

8. An apparatus for measuring and displaying the complex 
impedance of an object comprising: 
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frequency selection means for selecting at least one 
sinusoidal frequency; 

signal generating means for creating an excitation signal 
comprising a sinusoidal component at said selected sinusoidal 
frequency; 

response detection means for detecting a response signal 
arising in said object in response to application of said 
excitation signal; 

quadrature delay means for delaying said excitation 
signal by a predetermined time period to create a delayed 
excitation signal; 

first multiplier means for multiplying said excitation 
signal and response signal to form a first impedance signal; 

second multiplier means for multiplying said delayed 
excitation signal and response signal to form a second 
impedance signal; 

first averaging means for averaging said first impedance 
signal over time to form a first mean impedance signal; 

second averaging means for averaging said second 
impedance signal over time to form a second mean impedance 
signal; and 

display means for displaying said first and second mean 
impedance signals to represent said complex object impedance 
at said selected sinusoidal frequency. 

9, The apparatus described in Claim 8 further 
comprising: 

frequency memory means for storing a plurality of 
selected sinusoidal frequencies; and 

frequency sequencing means for stepping through said 
plurality of sinusoidal frequencies and sequentially 
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dispiaying said complex object impedance for each said 

5 sinusoidal frequency in turn. 

1C. The apparatus described in Claim 9 further 
2 comprising: 

calibration means for adjusting said first and second 
4 mean impedance signals in accordance with variations in 
excitation signal magnitude. 

11. The apparatus described in Claim 8 wherein: 

2 said sinusoidal frequency generating means , said 

response detection means, said delay means and said first and 

4 second multiplier means comprise digital electronic logic 
means * 

12. The apparatus described in Claim 11 wherein: 

2 said display means comprises digital display means for 

displaying numerical values of said complex body impedance. 

13 . The apparatus described in Claim 11 further 
2 comprising: 

frequency memory means for storing a plurality of 

4 selected sinusoidal frequencies; and 

frequency sequencing means for stepping through said 

6 plurality of sinusoidal frequencies and sequentially 
displaying said complex object impedance for each said 

5 sinusoidal frequency in turn. 
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14 . The apparatus described in Claim 13 further 
2 comprising: 

generator memory means for storing a plurality of digital 
4 words representing said excitation signal; and 

digital-to-analog conversion means for converting said 
6 digital words to an analog excitation signal. 

15. The apparatus described in Claim 14 wherein: 

2 said response detection means comprises analog-to- 

digital conversion means for sampling said response signal to 
4 . form a series of response samples. 

16. The apparatus described in Claim 15 wherein: 

2 said excitation signal comprises a sinusoidal electrical 

current and said response signal comprises a sinusoidal 
4 electrical voltage. 

^ 17. The apparatus described in claim 16 wherein: 
2 said object comprises living biological tissue. 

18. _ The apparatus described in Claim 17 further 
2 comprising: 

current limiting means for limiting said excitation 
4 signal current to a value that will cause no harm when applied 
to said living biological tissue. 

19. An apparatus for measuring and displaying at least 
2 one equivalent circuit component arising from a theoretical 

model of living biological tissue comprising: 
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4 frequency selection means for selecting at least one 
sinusoidal frequency; 

6 signal generating means for creating an excitation signal 

comprising a sinusoidal component at said selected sinusoidal 

5 frequency; 

response detection means for detecting a response signal 
10 arising in said object in response to application of said 

excitation signal; 
12 quadrature delay means for delaying said excitation 

signal by a predetermined time period to create a delayed 
14 excitation signal; 

first multiplier means* for multiplying said excitation 
16 signal and response signal to form a first impedance signal; 

second multiplier means for multiplying said delayed 
18 excitation signal and response signal to form a second 

impedance signal; 
20 first averaging means for averaging said first impedance 

signal over time to form a first mean impedance signal; 
22 second averaging means for averaging said second 

impedance signal over time to form a second mean impedance 
24 signal; 

first, memory means for storing said first and second mean 
2 6 impedance signals to represent said living biological signal 

impedance at said selected sinusoidal frequency; 
28 processing means for obtaining a plurality of said first: 

and second mean impedance signals representing said living 
2 0 biological tissue impedance at a plurality of selected 

sinusoidal frequencies; 
2 2 network analyzer means for inverting said plurality of 

living biological tissue impedances rc a plurality cf 
34 equivalent circuit element values; and 
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display means for displaying said plurality of living 
3 6 biological tissue equivalent circuit value. 
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